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Abstract 

In this article we perform a computational study of Polyrakis algorithms presented in 
[12,13]. These algorithms are used for the determination of the vector sublattice and the 
minimal lattice-subspace generated by a finite set of positive vectors of W k . The study 
demonstrates that our findings can be very useful in the field of Economics, especially in 
completion by options of security markets and portfolio insurance. 
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1 Introduction 



This paper provides computational methods for the determination of the vec- 
tor sublattice and the minimal lattice-subspace of M. k generated by a set B = 
{x\,X2, ■■■,x n } °f linearly independent positive vectors of M. k . In order to reach 
our goal the study of a vector-valued function (3 is further involved, which is de- 
fined and studied by I. Polyrakis in H12I13H . In particular, in H12I13II . it is proved 
that if B = {xi,x 2 , x n } is a set of linearly independent positive vectors of M. k 
and f3(i) = ^0f^ for each % = 1, 2, k, where r(i) = (x x (i), x 2 (i), x n {i)), 
R((3) denotes the range of the function (3 and K is the convex hull of R((3) then it 
holds, 
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(i) There exists an n— dimensional minimal lattice-subspace of IR fc containing B, 

(ii) The vector sublattice generated by B is an m— dimensional subspace of R fc . 



where n is the number of vertices of K and m is the number of different values of 



In IfTOll , the vector sublattice generated by a set of linearly independent positive 
vectors of R fe was used in order to provide an answer to the challenging question 
of market completion by options, in a discrete-time framework. Also, in [fTTTl it 
is proved that if for a set of linearly independent positive vectors of M fc there ex- 
ists a minimal lattice-subspace containing this set, then one can use the produced 
subspace for the study of a cost-minimization problem that is also known as the 
minimum cost portfolio insurance problem. From the computational point of view, 
several methods are presented in H7I8I91 in order to check whether a finite collection 
of linearly independent positive vectors(resp. functions) of IR fc (resp. C[a, b]) forms 
a lattice-subspace or a vector sublattice. In addition, these methods have been used 
for the calculation of the minimum-cost portfolio insurance. The main goals of this 
article are outlined as follows: 



(1) Exploit the theory of positive bases and derive new computational methods for 
the calculation of the completion by options of a two-period security market, 

(2) Calculate the minimal lattice-subspace generated by a finite collection of lin- 
early independent positive vectors. 



The proposed numerical solution, based upon previous works of the authors, pro- 
vide the exact solution to the problem of calculating the market completion by 
options and the minimum-cost insured portfolio. Also, a closer look at Example 
[12] can tell us that a manual procedure in order to determine the completion of the 
given security market can easily be a prohibited task. In the economics models that 
we are working the input data, as for example the payoffs of the different securities, 
are standard for a long time period and so, by our methods we determine the exact 
solution for these problems. 



In the present work, the numerical tasks have been performed using the high-end 
Matlab programming language. Specifically, the Matlab 7.4 (R2007b) environment 
J5H was used on an Intel(R) Pentium(R) Dual CPU T2310 @ 1.46 GHz 1.47 GHz 
32-bit system with 2 GB of RAM memory running on the Windows Vista Home 
Premium Operating System. 
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2 Lattice-subspaces and vector sublattices of R fc 

2.1 Preliminaries 

We first recall some definitions and notation from the vector lattice theory. Let 

R k = {x = (x(l), x(2), x(k))\x(i) G R, for each i}, where we view M fc as an 
ordered space. The pointwise order relation in M fc is defined by 

x < y if and only ifx{i) < y(i), for each i — 1, k. 

The positive cone of M. k is defined by = {x E M. k \x(i) > 0, for each i} and 
if we suppose that X is a vector subspace of IR fc then X ordered by the pointwise 
ordering is an ordered subspace of M fe with positive cone X + — X n R+. A point 
x G R fc is an upper bound (resp. lower bound) of a subset S C R fc if and only if 
y < x(resp. x < y), for all y E S. For a two-point set S = {x, y}, we denote by 
x V y(resp. x Ay) the supremum of 5 i.e., its least upper bound(resp. the infimum 
of 5 i.e., its greatest lower bound). Thus, x V y(resp. x A y) is the componentwise 
maximum(resp. minimum) of x and y defined by 

(x V y){i) = max{x(i), y(i)}((x A y){i) = mm{x(i), y(i)}), for alii = 1, ...,k. 

An ordered subspace X of R k is a lattice-subspace of R fe if it is a vector lattice 
in the induced ordering, i.e., for any two vectors x, y G X the supremum and the 
infimum of {x, y} both exist in X. Note that the supremum and the infimum of 
the set {x, y} are, in general, different in the subspace from the supremum and the 
infimum of this set in the initial space. An ordered subspace Z of R fc is a vector 
sublattice or a Riesz subspace of R k if for any x, y G Z the supremum and the 
infimum of the set {x, y} in IR fc belong to Z. 

Assume that X is an ordered subspace of IR fc and B = {bi,b 2 , b n } is a basis 
for X. Then B is a positive basis of X if for each x G X it holds that x is positive 
if and only if its coefficients in the basis B are positive. In other words, B is a 
positive basis of X if the positive cone X + of X has the form, 

n 

X + = {x = Aj6j|Aj > 0, for each i}. 

i=l 

Then, for any x = J27=i and y = Yh=\ A*A we nave x < ?/ if and only if 
Ai < f^i for each ? = 1, 2, n. 

Each element bi of the positive basis of X is an extremal point of X + thus a 
positive basis of X is unique in the sense of positive multiples. Recall that a nonzero 
element xq of X + is an extremal point of X + if, for any x E X,0 < x < x implies 
x = Xx for a real number A. The existence of positive bases is not always ensured, 
but in the case where X is a vector sublattice of R k then X has always a positive 
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basis. Moreover, it holds that an ordered subspace of M. k has a positive basis if and 
only if it is a lattice-subspace of IR fc . If B = {bi, b 2 , b n } is a positive basis for a 
lattice-subspace (or a vector sublattice) X then the lattice operations in X, namely 
xVy for the supremum and x A y for the infimum of the set {x, y} in X, are given 
by 

n n 

X ^V = H max {^ ; and x A y = ^ min{Ai, fii}b h 

i=l i=l 

for each x = I]™ =1 Ajfoj, y = J27=i fak E X. A vector sublattice is always a lattice- 
subspace, but the converse is not true as shown in the next example, 

Example 1 Let X = [x 1 ,x 2 , x s ] be the subspace of IR 4 generated by the vectors 
x\ = (6, 0, 0, 1), x 2 = (6, 4, 0, 0), x 3 = (8, 4, 2, 0). An easy argument shows that 
the set B = 62, 63} where 



61 




2 


2 


b 2 




12 


8 






6 


1 



forms a positive basis of X therefore X is a lattice-subspace of IR 4 . On the other 
hand, let us consider the vectors y\ = 2xi + x 2 = (18, 4, 0, 2) and y 2 = x 3 — x 2 = 
(2, 0, 2, 0) of X. Then, y x V y 2 = (18, 4, 2, 2) and since y a = \b 2 + 2fo 3 , y 2 = b x , it 
follows that y\Vy 2 = b± + |6 2 + 26 3 = (20, 4, 4, 2). From the definition of a vector 
sublattice we have that X is a vector sublattice of IR 4 if for each x, y E X it holds 
x\/y = xVy G X and x Ay = x A y E X . Therefore, since there are two elements 
y±,y 2 E X such that y 1 V y 2 7^ yi Vy 2 it follows that X is not a vector sublattice of 
IR 4 . 

For an extensive presentation of lattice-subspaces, positive bases and vector sub- 
lattices, we refer to H1I11I12I13II . for computational methods in positive bases the- 
ory we refer to H7I8I9B and for several applications in the theory of finance the reader 
may refer to I2I3I4I101 . 



2.2 The mathematical problem 

In this section we shall present the results of H12I13H for the construction of a min- 
imal lattice-subspace of lR fe containing a linearly independent subset of Ri together 
with the construction of the vector sublattice of IR fc that this subset generates. In par- 
ticular, given a collection of linearly independent, positive vectors, Xi, x 2 , x n of 
]R fc then the basic tool for our study is a special function, named basic function, that 
this collection of vectors defines. This function was first introduced and studied in 
|fi~2|. In our analysis we shall use the notation introduced in [12], so let us denote 
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by r the function r : {1 



2, k} -> M fe such that 



r(i) = (xi(i),x 2 («)j ■••> a; n(0) 



and by /3 the function (3 : {1, 2, k} -)■ IR fc such that 



for each i £ {1,2, A;} with ||r(i) ||i ^ 0. We shall refer to (3 as the basic function 
of the vectors x\, x 2 , x n . The set 



is the range of the basic function and the cardinal number, cardR((3), of R{0) is 
the number of different elements of R(/3). Let cardR((3) = m then n < m < k 
and by K we shall denote the convex hull of R(j3) which is, as the convex hull of 
a finite subset of R fc , a polytope with d vertices and each vertex of K belongs to 
R(/3). It is clear that n < d < m. 

An essential issue for our analysis are the conditions under which a collection of 
linearly independent, positive vectors x\, x 2l x n of IR fc can be used to derive a 
minimal lattice-subspace and a vector sublattice containing these vectors. In IfTBl , 
theorem 3.19 is a criterion for lattice-subspaces and vector sublattices and provides 
a full answer on the topic on the basis of describing the geometry of this problem. 
In the case where X = [xx, x 2 , x n ] is a lattice-subspace or a vector sublattice of 
]R fc then the theorem determines a positive basis in X while, in the opposite case, 
theorem 3. 19 provides a minimal lattice-subspace and a vector sublattice containing 
X. In order to state theorem 3.19, let us consider R{(3) = {Pi, P 2 , P m } such that 
the first n vertices Pi, P 2 , P n are linearly independent and Pi, P 2 , P^ are the 
vertices of K, n < d < m. Also, A T denotes the transpose matrix of a matrix A. 

Theorem 2 [13, Theorem 3.19]. Suppose that the above assumptions are satisfied. 
Then, 

(i) X is a vector sublattice ofM k if and only if R{j3) has exactly n points (i.e., 
m = n). Then a positive basis b\,b 2 , b n for X is defined by the formula 



where A is the n x n matrix whose ith column is the vector P, for each 
i — 1, 2, m. 

(ii) X is a lattice-subspace ofM. k if and only if the polytope K has n vertices (i.e., 
d = n). Then a positive basis &i, b 2 , . . . , b n for X is defined by the formula 



R((3) = {(3(i)\i = 1,2,.., A;, with ||r(i)||i ^ 0} 



Oi, b 2 , .., b n ) T = A 1 {x 1 ,x 2 , .., x n ) T 



(61,62 




■^1 , %2 j • • • 5 ■En) i 
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where A is the n x n matrix whose ith column is the vector Pi, for each 
i = 1,2,..., d. 
(Hi) Let m > n. Ifl s = j3~ l (P s ), and 

x s = Yl H r 00lli e i' s = n + l,n + 2, . . . ,m, 

then 

Z \X\ , X 2 j • • • ; X n j -En+l j • • • ) •Em\ 

is the vector sublattice generated by X\ , x 2 , . . . , x n and dim Z = m. 
(iv) Let d > n.If^: D(f3) -»■ R + , i = 1,2,..., d such that Yh=i = 1 an d 
= Y,i=i€i{j)Pifor each j e D(/3), and x n+i ,i = 1,2, ... ,d - n, are 
the following vectors ofR k : 

Xn+i= ^n+i(j)\\r(j)\\iej, 

then 

Y \p^l j • • • ; %n ■> 1 ) ■ ■ ■ ) ■Ed\ 

is a minimal lattice-subspace ofR k containing X\, x 2 , ■ ■ ■ , x n and dim Y — d. 

In (H, a new numerical package is provided so that, for a variety of dimensions 
and subspaces, we are able to check conditions (i), (ii) of Theorem [2l In addition, 
in the case where X = [x±, x 2 , ■ ■ ■ , x n ] is a vector sublattice or a lattice-subspace 
the proposed method, in [9|, provides the positive basis of X therefore a complete 
description of the vector sublattice or the lattice-subspace is given. 

In the case where X is not a vector sublattice or a lattice-subspace, an essential 
issue for our analysis is the construction of a powerful and efficient package in 
order to calculate the vector sublattice and a minimal lattice-subspace containing 
X by using conditions (Hi), (iv) of Theorem|2l 



2.3 The algorithm 



(1) Determine the function /3 as well as the range, R(f3), of (3. 

(2) Compute the number m = cardR((3), and the number d of vertices of the 
polytope K. 

(3) If n — m (vector sublattice case) or n = d (lattice-subspace case) then, deter- 
mine a positive basis of X. 

(A) If m > n, then Z = [xi,x 2 , . . . , x n , x n+ i, . . . , x m ] is the vector sublattice 

generated by x u x 2 , ...,x n , where x s = J2 i€ i s IkWIlie* and I s = /? _1 (P S ) 

for each s = n + l,n + 2, . . . ,m . 
(5) If d > n, then Y — [xi,...,x n , x n+ \, . . . , Xd] is a minimal lattice-subspace of 

R k containing xi, x 2 , . . . , x n , where the vectors x n+i = J2jeD(p) £n+i(j) \\r(j) || \e h 

% = 1, . . . , d — n, were defined in (iv) of Theorem[2l 
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Note that, the steps of this algorithm are based upon Theorem[2] In the following 
section, we present the translation followed by the implementation of this algorithm 
in M fc within two Matlab-based functions named SUB 1 at and MINI at . These 
functions, as we shall see in Subsection 13. 1[ provides an important tool in order to 
investigate minimal lattice-subspaces and vector sublattices of IR fc generated by a 
set of linearly independent, positive vectors Xi, x%, x n of R. k . 



3 The computational method 

3.1 Method presentation and examples 

Our proposed numerical method is based on the introduction of two functions, 
namely SUB 1 at and MINI at (see Appendix), that enable us to perform fast test- 
ing for a variety of dimensions and subspaces. Also, both of these functions are 
using the function SUB 1 at SUB from [9], in order to calculate a positive basis. For 
the sake of completeness, we provide the SUBlatSUB function in the Appendix. 
Recall that, the numbers n,m,d,k denote the dimension of X, the cardinality of 
R(/3), the number of vertices of the convex hull of R(f3) and the dimension of the 
initial Euclidean space, respectively. 

The function SUBlat first determines R((3) and the number m = cardR((3) and 
then the n linearly independent vertices Pj, % — 1, . . . , n of the poly tope K. Finally, 
the program calculates the subsets I s = f3~ 1 (P s ) and then the vectors 

x s = J2 ll r Wlli e ;> 

for each s — n + l,n + 2, . . . , m. Therefore, since by theorem [2] (in) the subspace 

Z [Xl, 3^2, ••• j X n: X n _|_i, . . . , X m ] 

is the vector sublattice generated by the given vectors x\, x 2 , ■ ■ ■ , x n the program 
responds, by using the function SUBlatSUB of O, with two mx k matrices. The 
rows of the first matrix are the vectors 

while the rows of the second matrix are the vectors bi, . . . , b m of the positive basis 
for the subspace Z. 

In order to determine a minimal lattice-subspace generated by the given collec- 
tion of vectors we use the function MINI at . The MINI at function first deter- 
mines R((3) and then the vertices of the polytope K. The correct performance of 
the MINI at function requires the use of the convhulln Matlab function which is 
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based on Qhull (For information about Qhull see http://www.qhull.org/). In order to 
determine a minimal lattice-subspace of R. k containing x\, x 2 , ■ ■ ■ , x n the program 
calculates the vectors £j : D((3) — > R + , i = 1, 2, . . . , d by solving underdeter- 
mined nxrf linear systems subject to the inequalities & > ,i = 1, . . . n. Finally, 
the program defines the vectors, 



x n+i = £n+iU)\\ r U)h e j> i = l,2,...,d-n, 

jeD(/3) 



and since the subspace 



is a minimal lattice-subspace containing a?i, x 2 , . . . , x n the program provides, by 
using the function SUBlatSUB of flU, two d x k matrices. The rows of the first 
matrix are the vectors 



while the rows of the second matrix are the vectors bi, b 2 , ■ ■ . , b d of the positive 
basis for the subspace Y. 



Consequently, let xi, x 2 , x n be a collection of linearly independent, positive 
vectors of M fc , then we construct a matrix 5 whose columns are the vectors of the 
given collection and then we apply the functions SUBlat and MINlat on that 
matrix as follows, 



[ Vector Sub lattice , Posit iveba sis ] = SUBlat (B) 
[Minimal lattice sub space, Posit iveba s is ] =MINlat (B) 



In order to illustrate the most important features of SUBlat and MINlat , we 
reproduce two examples featured in H7I13H . 
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Example 3 Consider the following 10 vectors xi, x%, xio in IR 17 , 
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In it is provided that the subspace X = [xi,x%,. . . ,xio) is not a lattice- 
subspace of M 17 . In order to calculate the vector sublattice that X generates as 
well as a minimal lattice-subspace containing X we invoke the SUBlat and the 
MINI at functions by typing in the command window of the Matlab environment: 

>> [ VectorSublatt ice, Posit ivebasis] =SUBlat (B) 

>> [Minimal lattice sub space, Posit ivebasis ] =MINlat (B) 

where B denotes a matrix whose columns are the vectors x±, x 2 , ■ ■ ■ , x w . The re- 
sults, then, are as follows: 



VectorSublattice = 
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Minimallatticesubspace = 
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Example 4 Consider the following 4 vectors x 1 ,x 2 ,x 3 , x 4 in M 7 , 



Xi 




12 10 114 


x 2 




11110 2 


x 3 




2 10 1112 


X4 




10 1110 



where following the same procedure, as before, one gets 



VectorSublattice = 

12 10 114 
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Remark 5 In the beginning of the present section, we mentioned that in order to 
determine the vectors & : D(f3) — > IR + , i = 1,2, ... ,d the function MINlat re- 
sponds by solving k underdetermined nxd linear systems subject to the inequalities 
6 > ,i = 1, . . . , n. Since, the solution to such systems is not unique then, if p 
denotes a solution corresponding to one of these systems, then y = p + h, where h 
is an arbitrary vector from the null space, is a solution too. Thus, a minimal lattice- 
subspace containing the vectors Xi, x 2 , x n is not unique. In |fT3ll . Example 3.21, 
it is proved that for the collection of vectors of Example |4] there exist two minimal 
lattice-subspaces Y, Y' such that 

• Y^Y', 

• the space Y D Y' is not a lattice-subspace, 

• and Y, Y' are not subspaces of Z (Z denotes the generated vector sublattice). 



3.2 Comparison Results 



For the purpose of monitoring the performance, in this section we present the 
execution times of the proposed methods (SUB 1 at, MINlat ) for various col- 
lections of vectors and dimensions. For this purpose we have used the function 
testMINlat , testSUBlat (see Appendix) in order to test 50 full rank matri- 
ces for each rank n, n — 3, 30. The cumulative results are presented in Table 1 
while, the time responses have been recorded using the Matlab function profile. 

It is evident, from Table 1, that the proposed numerical methods, based on the 
introduction of the SUB 1 at and the MINlat function, enable us to perform fast 
estimations for a variety of dimensions. 
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Table 1 



Results for 50 tested full rank matrices for each rank 


n, n = 3, 30. 




Rank 


SUBlat 


MINlat 


Rank 


SUBlat 


MINlat 




(Total time in seconds) 


(Total time in seconds) 




(Total time in seconds) 


(Total time in seconds) 


3 


0.553 


0.689 


17 


2.812 


1.528 


4 


0.569 


0.705 


18 


2.994 


1.732 


5 


0.584 


0.742 


19 


3.333 


1.885 


6 


0.696 


0.822 


20 


3.594 


2.222 


7 


0.856 


0.831 


21 


4.084 


2.428 


8 


0.973 


0.890 


22 


4.219 


2.773 


9 


1.141 


0.903 


23 


4.569 


3.175 


10 


1.307 


1.012 


24 


4.817 


3.695 


11 


1.642 


1.047 


25 


5.149 


4.306 


12 


1.667 


1.108 


26 


5.712 


4.815 


13 


1.890 


1.151 


27 


5.910 


5.741 


14 


2.059 


1.185 


28 


6.229 


6.566 


15 


2.314 


1.270 


29 


6.675 


7.493 


16 


2.765 


1.365 


30 


7.003 


9.226 



4 Applications in economies with incomplete asset markets 



The theory of vector sublattices and lattice-subspaces has been extensively used 
in the last years in Mathematical Economics, especially in the areas of incomplete 
markets and portfolio insurance. Recently, in IfTOl , a new approach to the problem 
of completion by options of a two-period security market was introduced, which 
used the theory of positive bases in vector sublattices of M. k (as described in Ifl3l0 . 
Also, if xi, x%, . . . ,x n denotes a collection of linearly independent, positive vec- 
tors, of IR fc then if X = [xx, x 2 , ■ ■ ■ , x n ) is a lattice subspace or, in the contrary case, 
if there exists a minimal lattice-subspace containing X then a solution to a cost 
minimization problem known as minimum-cost portfolio insurance always exists. 
In this section we shall discuss these interconnections and briefly describe the the- 
oretical background on this subject. A full theoretical development is available in 
H2I3I4I7I8I10I1 1I12I13H . We will be also presenting how one can use the proposed 
functions, SUBlat and MINlat, in order to determine the completion of security 
markets as well as the solution of the minimum-cost portfolio insurance problem. 



4.1 Completion of security markets 



The problem of completion by options of a two-period security market in which 
the space of marketed securities is a subspace of R fc has been studied in IfTOl . The 
present study involves vector sublattices generated by a subset B of IR fc of positive, 
linearly independent vectors, so we shall provide a computational solution to this 
problem by using the SUBlatfunction in order to provide the generated vector 
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sublattice Y as well as a positive basis for Y. 

Let us assume that in the beginning of a time period there are n securities traded 
in a market. Let S = {1, .... k} denote a finite set of states and Xj G IR+ be the pay- 
off vector of security j in k states. The payoffs Xi, x 2 , x n are assumed linearly 
independent so that there are no redundant securities. If 9 = [Ox, 2 , 6 n ) G IR n is 
a non-zero portfolio then its payoff is the vector 

n 

T(9) = y £e i x i . 

i=X 

The set of payoffs of all portfolios is referred as the space of marketed securities 
and it is the linear span of the payoffs vectors x\, x 2 , x n in M fe which we shall 
denote it by X, i.e., 

X = [xx,x 2 , ...x n }. 

For any x, u G M fc and any real number a the vector c u (x, a) = (x — au) + is the 
call option &rx&p u (x, a) = {au — x) + is the put option of x with respect to the strike 
vector u and exercise price a. 

Following the terminology of IfTOl , let U be a fixed subspace of M. k which is called 
strike subspace and the elements of U are the strike vectors. Then, the completion 
by options of the subspace X with respect to U is the space Fu{X) which is defined 
inductively as follows: 

• Xx is the subspace of IR fc generated by 0\, where Ox = {c u (x, a)\x G X,u G 
U, a G M.}, denotes the set of call options written on the elements of X, 

• X n is the subspace of M fc generated by O n , where O n = {c u {x,a)\x G X n _i,w G 
U, a G M.}, denotes the set of call options written on the elements of X n _ l5 

. Fu(X) = U™ =1 X n . 

The completion by options F V {X) of X with respect to U is the vector sublattice 
of IR fc generated by the subspace Y — X U U. The details are presented in the next 
theorem, 

Theorem 6 / fTQl Theorem 3] In the above notation, we have 

(i) ^ C Xx, 

{ii) Fjj{X) is the sublattice S{Y) ofM k generated by Y, and 
(Hi) ifU CI, then Fjj(X) is the sublattice ofW k generated by X. 

Definition 7 / f70l Definition 10] Any set {yx, y 2 , ■ ■ ■ ,y r } of linearly independent 
positive vectors ofM h such that Fu(X) is the sublattice of~R k generated by {yx, y 2 , . . . , 
is a basic set of the market. 

Theorem 8 [10. Theorem 11] Any maximal subset {yx, y 2 , . . . ,y r } of linearly in- 
dependent vectors of A is a basic set of the market, where A = {xf , x~[ , . . . , x+ , x~ }, 



13 



ifU C X and A = {x± , x x , . . . , x+, x n , uf, u x , . . . , u£, u d }, ifU C X 

Definition 9 £70] Definition 12] The space of marketed securities X is complete 
by options with respect to U if X = Fu(X). 

From theorem [2] and definition [9] it follows, 

Theorem 10 / liOl Theorem 13] The space X of marketed securities is complete by 
options with respect to U if and only ifU^X and cardR{j3) = n. 

Theorem 11 / fTQl Theorem 14] The dimension of Fjj(X) is equal to the cardinal 
number of R(/3). Therefore, Fu(X) = M fc if and only if cardR(j3) = k. 

In order to apply our method to the problem of completion of security markets 
we present the following example which is formerly featured in ifTOl . 

Example 12 / fTQl Example 16] Suppose that in a security market, the payoff space 
is R 12 and the primitive securities are: 

x x = (1,2,2,-1,1,-2,-1,-3,0,0,0,0) 
x 2 = (0,2,0,0,1,2,0,3,-1,-1,-1,-2) 
x 3 = (1,2, 2, 0,1, 0,0, 0,-1, -1,-1, -2) 
and that the strike subspace is the vector subspace U generated by the vector 

u=(l, 2,2, 1,1,2, 1,3, -1,-1, -1,-2). 

Then, a maximal subset of linearly independent vectors of {xf , x± , x\ , x 2 , x% , x% , uf , u j~ } 
can be calculated by using the following codeQ: 

>>XX = [max (X, zeros (size (X) ) ) ; max (-X, zeros (size (X) ) ) ] ; 
>>S = rref (XX' ) ; 
>> [I, J] = f ind (S) ; 

>>Linearindep = accumarray (I, J, [rank (XX) , 1] , @min) ' ; 
>>W = XX (Linearindep, :) ; 

where X denotes a matrix whose rows are the vectors xi, x 2 , x 3l u. We can de- 
termine the completion by options of X i.e., the space Fjj(X), with the SUBlat 
function by using the following code: 

>> [ VectorSublatt ice, Posit ivebasis] =SUBlat (W ) 
The results then are as follows 



The rref function is a Matlab function that produces the reduced row echelon form of a 
given matrix by using Gauss Jordan elimination with partial pivoting (cf. 1151"). 
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VectorSublattice = 

122010000000 
020012030000 
122112130000 
000000001112 
204000000000 

Positivebasis = 

000000001112 
000100100000 
000004060000 
408000000000 
060030000000 



4.2 Minimum cost portfolio insurance 



In this section we shall, briefly discuss an investment strategy called minimum- 
cost portfolio insurance as a solution of a cost minimization problem. We will be 
also presenting how one can use the previous results for the calculation of the min- 
imal lattice-subspace, in order to calculate the minimum-cost insured portfolio. In 
our model we use a method of comparing portfolios called portfolio dominance 
ordering. This ordering compares portfolios by means of the ordering of their pay- 
offs. Under this consideration we are able to use the order structure of the payoff 
space together with the theory of lattice-subspaces. In what follows we shall use 
the notation introduced in IfTTTl . 

The model of security markets we study here is extended over two periods, 
namely period and period 1. We assume n securities labeled by the natural num- 
bers 1, 2, n, acquired during the period and that these n securities are described 
by their payoffs at date 1. The payoff of the ith security is in general a positive ele- 
ment xi of an ordered space E which is called payoff space. In addition, we assume 
that the payoffs X\ , x 2 , ■ ■■ , x n are linearly independent so that there are no redundant 
securities and that the securities have limited liability which ensures the positivity 
of xi, X2, x n . In [fTTI . it is assumed that E is the space of real valued continuous 
functions C(Vt) defined in a compact, Hausdorff topological space Vt. A portfolio 
is a vector 9 = (9 1 , 9 2 , 9 n ) of M n where 9i is the number of shares of the ith se- 
curity. The space W 1 is then known as portfolio space. If 8 = 6 2 , n ) e M n 
is a non-zero portfolio then its payoff is the vector 

R{e) = y £e i x i ec{n). 

i=i 

The operator R is one-to-one and is called the payoff operator. The pointwise or- 
dering in C(Q), induces the partial ordering > R in the portfolio space W 1 and is 
defined as follows: For each 9, 4> G W 1 we have 

>r <P, if and only if R{9) > R{4>). 
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This ordering is known as the portfolio dominance ordering. The set of payoffs of 
all portfolios, or the range space of the payoff operator, is the linear span of the 
payoffs vectors xi,x 2 , x n in C(O) which we shall denote it by Ai, i.e., 

M = [x 1} x 2 , ...,x n ]. 

The subspace M. of C(Vt) is called the asset span of securities or the space of 
marketed securities. 

Let us assume that p = (pi,p 2 , ...,p n ) £ M fc is a vector of security prices and 
are two portfolios. Then, the insured payoff on the portfolio 8 = (8 1 ,8 2 , 8 n ) at 
the "floor" </> and in the price p is the contingent claim R(8) V R{4>). 

The solution of the following cost minimization problem is referred to as the 
minimum-cost insured portfolio, or a minimum-cost insurance of the portfolio 8 at 
the floor (f) and in the price p, 

min p ■ r] 

T)GK fc 

subject to 

R{rf) > R(8) V R((j>). 

In [fTTI it is proved that if the payoff space is contained in a minimal lattice- 
subspace of C(f2) then a minimum-cost insurance of the portfolio 8 always exists. 
The details are included in the next theorem: 

Theorem 13 / [771 Theorem 15] If the payoff space X is contained in a finite- 
dimensional minimal lattice-subspace Y of C(Q) and the sum of the payoff vec- 
tors Xi is strictly positive, then a minimum-cost insurance of the portfolio 8 at the 
floor cf) and in the price p, exists and it is determined by solving the corresponding 
minimization problem. 

In order to apply our method, MINI at, we consider that = {1, 2, n}. Then, 
it is evident that C(fi) = M n . Therefore, in view of Theorem [T3l and by using 
the MINI at function we are able to determine the minimum-cost insurance of the 
given portfolio 8 at the floor cf) and in the price p. 



5 Conclusions 

In this paper, new computational methods in order to determine vector sublattices 
and minimal lattice-subspaces of M. k are presented. In order to reach our goal the 
study of a vector- valued function (3 is further involved by introducing two Matlab 
functions, namely SUB 1 at, and MINI at. The results of this work can give us an 
important tool in order to study the interesting problems of completion by options 
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of a two-period security market in which the space of marketed securities is a sub- 
space of IR fc and in portfolio insurance. The experiment results, in subsection 13 .2[ 
show that our algorithm performs well. Also, note that the algorithm |2.3| determines 
the exact solution to the problems of completion by options and the minimum-cost 
insured portfolio. Finally we are convinced that, from the mathematical point of 
view, the proposed algorithm can be further analyzed independently, in terms of 
formal numerical analysis. 



6 Appendix 

The SUB 1 at function 

function [ Sublattice, Posit ivebasis ] = SUBlat (B) 
%SUBlat (B) provides the vector sublattice generated by 
%a given finite collection of positive, linearly 
%independent vectors of R"n 

%B denotes the matrix whose columns are the given vectors 
[N, M] = size (B) ; 
Id = eye (N) ; 
for i = 1:N, 

if norm (B (i, : ) , 1) ~=0, 
Test(i, :) = 1/ norm (B(i,:),l)*B(i,:); 

end 
end 

Matrix = Test; 

[BB,m,n] = unique (Matrix, ' rows ') ; 
Index = 1:N; 
S = rref (BB' ) ; 

[I, J] = f ind (S) ; 

Linear indep = accumarray ( I , J, [ rank (BB) , 1 ] , @min) ' ; 
mm = length (m) ; 
nn = length (n) ; 
Indexl = 1 : mm; 

Index2 = setdiff ( Indexl , Linearindep) ; 
YY = sum (B, 2) ' ; 

TTT = setdiff ( Index, m (Linearindep) ) ; 

KK = Id (TTT, : ) ; 

TT = YY (1, TTT) ' ; 

T = diag (TT) *KK; 

K = zeros (N) ; 

K (TTT, : ) = T; 

Vec = zeros (mm-M, N) ; 

if mm < nn, 
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for i = 1 : length ( Index2 ) , 

DD = strmatch ( Index2 ( i ), n, ' exact ') ' ; 
R = length (DD) ; 
if R >= 2, 

Vector = sum (K (DD, : ) ) ; 
else 

Vector = K (DD, : ) ; 

end 

Vec(i,:) = Vector; 

end 

[a,b] = find (Vec) ; 
Vectors = Vec (unique (a) ,:) ; 
Sublattice = [B' ; Vectors] ; 
Positivebasis = SUBlatSUB (Sublattice' ) ; 

else 

KKK = unique (K, ' rows ') ; 
[II, JJ] = find (KKK) ; 
Vectors = KKK (unique ( I I ) , : ) ; 
Sublattice = [B' ; Vectors] ; 
Positivebasis = SUBlatSUB (Sublattice' ) ; 

end 

The MINI at function 

function [MLS , Posit ivebasis ] = MINlat (B) 
%MINlat (B) provides the minimal latt ice-subspace 
%generated by a given finite collection of positive, 
%linearly independent vectors of R"n 

%B denotes the matrix whose columns are the given vectors 
[N,M] = size (B) ; 
for i=l : N, 

if norm (B (i, : ) , 1) ~=0, 

Test(i, :) = 1 /norm (B ( i , : ) , 1 ) *B ( i , : ) ; 

end 
end 

Matrix = Test; 
[ii,jj] = find (Matrix) ; 
Matrixl = Matrix (unique (ii) , : ) ; 
BB = unique (Matrixl ,' rows ') ; 
Ml = rank (bsxfun (@minus, BB, BB (1, : ) ) ) ; 
if MKM, 

Utrans = bsxf un ( @minus , BB, BB ( 1 , : ) ) ; 
Rot = orth (Utrans ') ; 
Uproj = Utrans*Rot; 
Tri = convhulln (Uproj ) ; 
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VIndex = unique (Tri ( : ) ) ; 
P = BB (VIndex, : ) ' ; 
Q = length (VIndex) ; 
else 

VIndex = unique (convhulln (BB) ) ; 
P = BB (VIndex, : ) ' ; 
Q = length (VIndex) ; 
end 

Test = zeros (N, M) ; 
for i=l : N, 

Test(i, :) = 1 /norm (B ( i , : ) , 1 ) *B ( i , : ) ; 

end 

Matrix = Test; 
BBB=Matrix' ; 
R = Q-M; 

Sol = zeros (Q, R) ; 
for i = 1:N, 

Sol(:,i) = lsqnonneg (P, BBB ( : , i) ) ; 

end 

Solutions = Sol; 
Norms = sum(B,2)'; 
Testl = zeros (R, N) ; 
for i = 1:R, 
Indexl = M+i; 

D = Solutions ( Indexl ,:). *Norms*eye (N) ; 
Testl (i, : ) = D; 
end 

Minlatsub = [B'; Testl]; 
MLS = Minlatsub; 

Positivebasis = SUBlat SUB (MLS ' ) ; 

The SUBlat SUB function 

function [posit ivebasis , dimensions ] = SUBlatSUB (A) 

%SUBlatSUB (A) provides the vector sublattice 

%or the latt ice-subspace of a given finite collection 

%of positive, linearly independent vectors of R"n 

%A denotes the matrix whose columns are the given vectors 

if any (any (A<0) ) ~=0, 

error ('the initial matrix must have positive elements') 

end 

[N,M] = size (A) ; 
if rank (A) ~=M, 

error ('the given vectors are linearly dependent') 

end 
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for i=l : N, 
if norm (A (i, : ) , 1) ~=0, 
Test(i,:) = 1 /norm (A ( i , : ) , 1 ) *A ( i , : ) ; 

end 
end 

matrix = Test; 
[ii,jj] = find (matrix) ; 

matrixl = matrix (unique (ii) , : ) ; 

u = unique (matrixl rows ') ; 

m = length (u ( : , 1 ) ) ; 
if M == m, 

disp ('vector sublattice' ) 

posit ivebasis = inv(u')*A'; 

dimensions = [M m N] ' ; 
else 

ml = rank (bsxf un ( @minus , u, u ( 1 , : ) ) ) ; 
if ml<M, 

utrans = bsxfun (@minus,u,u (1, : ) ) ; 

rot = orth (utrans ') ; 

uproj = utrans*rot; 

tri = convhulln (uproj ) ; 

d = length (unique (tri (:))); 

if d == M, 

basis = inv (u (unique (tri (:)),:)') *A' ; 

disp (' lattice-subspace' ) 

posit ivebasis = basis; 

dimensions = [M m d N] ' ; 

else 

disp ('not a lattice-subspace') 
dimensions = [M m d N] ' ; 
posit ivebasis= [ ] ; 
end 
end 
end 

The test SUB lat function 

function testSUBlat = test SUBlat (k, j ) 

%k is the the dimension of the Euclidean space R"k 

%j is the number of the tested matrices 

for i = 1 : j , 

A = rand (k+2, k) ; 

testSUBlat = SUBlat (A) ; 

end 
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The testMINlat function 

function testMINlat = testMINlat (k, j ) 

%k is the the dimension of the Euclidean space R"k 

%j is the number of the tested matrices 

for i = 1 : j , 

A = rand (k+2, k) ; 

testMINlat = MINlat (A) ; 

end 
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